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ABSTRACT 

The observed number counts of quasars may be explained either by long-Uved activity within rare 
massive hosts, or by short-hved activity within smaller, more common hosts. It has been argued that 
quasar lifetimes may therefore be inferred from their clustering length, which determines the typical mass 
of the quasar host. Here we point out that the relationship between the mass of the black-hole and the 
circular velocity of its host dark-matter halo is more fundamental to the determination of the clustering 
length. In particular, the clustering length observed in the 2dF quasar redshift survey is consistent with 
the galactic halo - black-hole relation observed in local galaxies, provided that quasars shine at ^ 10- 
100% of their Eddington luminosity. The slow evolution of the clustering length with redshift inferred in 
the 2dF quasar survey favors a black-hole mass whose redshift-independent scaling is with halo circular 
velocity, rather than halo mass. These results are independent from observations of the number counts 
of bright quasars which may be used to determine the quasar lifetime and its dependence on redshift. 
We show that if quasar activity results from galaxy mergers, then the number counts of quasars imply an 
episodic quasar lifetime that is set by the dynamical time of the host galaxy rather than by the Salpeter 
time. Our results imply that as the redshift increases, the central black-holes comprise a larger fraction 
of their host galaxy mass and the quasar lifetime gets shorter. 

Subject headings: cosmology: theory - cosmology: observations - quasars: general 



1. INTRODUCTION 

The Sloan Digital Sky Survey (York et al. 2000) and 
the 2dF quasar redshift survey (Groom et al. 2001a) have 
measured redshifts for large samples of quasars, and de- 
termined their luminosity function over a wide range of 
redshifts (Boyle et al. 2000; Fan et al. 2001a,b; Fan et 
al. 2003). The 2dF survey has also been used to constrain 
the clustering properties of quasars (Groom et al. 2001b). 
It has been suggested that quasars have clustering statis- 
tics similar to optically selected galaxies in the local uni- 
verse, with a clustering length i?o ~ 8Mpc. The large 
sample size of the 2dF quasars also provided clues about 
the variation of clustering length with redshift (Groom et 
al. 2001b) and apparent magnitude (Groom et al. 2002). 
Quasars appear more clustered at high redshift, although 
with a relatively mild trend. There is also evidence that 
more luminous quasars may be more highly clustered. 

As pointed out by Martini & Weinberg (2001) and 
Haiman & Hui (2001), the quasar correlation length de- 
termines the typical mass of the dark matter halo in 
which the quasar resides. One may therefore derive the 
quasar duty-cycle by comparing the number density of 
quasars with the density of host dark matter halos. The 
quasar lifetime then follows from the product of the duty- 
cycle and the lifetime of the dark-matter halo in between 
major mergers, although there is a degeneracy between 
the lifetime and the quasar occupation fraction or beam- 
ing. Preliminary results suggested quasar lifetimes of 

~ 10^ — 10^ years, consistent with the values determined 
by other methods (see Martini 2003 for a review), includ- 
ing the transverse proximity effect (Jakobsen et al. 2003) 
and counting arguments relative to the local population 



of remnant supermassive black holes (SMBHs; see Yu & 
Tremaine 2002). Kauffmann & Haehnelt (2002) used a de- 
tailed semi-analytic model (Kauffmann & Haehnelt 2000) 
to predict the correlation length of quasars and its evolu- 
tion with redshift. They found that their model reproduces 
the correct correlation length as well as its redshift evolu- 
tion for present-day quasar lifetimes of tq ~ 10*" years. 

In this work we argue that the correlation length of 
quasars is fundamentally determined by the relation be- 
tween the masses of SMBHs and their host galactic halo 
rather than by the quasar lifetime. We find that the am- 
plitude of the correlation length depends on the product 
of the SMBH-halo relation and the typical fraction of the 
Eddington luminosity at which quasars shine. Moreover, 
we show that the evolution of the correlation length is 
sensitive to how the SMBH - halo relation evolves with 
redshift, and therefore to the physics of SMBH formation 
and quasar evolution. 

The paper is organized as follows. In § 2 and § 3 we dis- 
cuss the calculation of the correlation function of quasars 
with a particular apparent magnitude B, and the differ- 
ent scenarios for the evolution of the SMBH ~ galaxy halo 
relation. In § 4 we compare fiducial model correlation 
functions to the results of the 2dF quasar redshift survey 
(Groom et al. 2001b, c). We find that the shallow evolution 
of the correlation length implies that SMBH mass has a 
redshift-independent scaling with the circular velocity of 
the host dark matter halo rather than with its mass. The 
ranges of the normalizations in the SMBH - galaxy halo 
relation and of the fraction of Eddington allowed by ob- 
servations of the observed quasar correlation function are 
explored in § 5. In § 6 we examine the quasar lifetime 



^ University of Melbourne, Parkville, Victoria, Australia 
Harvard-Smithsonian Center for Astrophysics, 60 Garden St., Cambridge, MA 02138 

1 



2 



by requiring that observational constraints from both the 
correlation and luminosity functions be satisfied simulta- 
neously. Finally, in § 8 we summarize our results and dis- 
cuss a very simple, physically motivated model that satis- 
fies all constraints with no free parameters. Throughout 
the paper we adopt the set of cosmological parameters de- 
termined by the Wilkinson Microwave Anisotropy Probe 
(WMAP, Spergel et al. 2003), namely mass density pa- 
rameters of fim = 0.27 in matter, fib = 0.044 in baryons, 
JIa = 0.73 in a cosmological constant, and a Hubble con- 
stant of Ho = 71 kms-i Mpc'^ 

2. THE CORRELATION FUNCTION OF QUASARS 

The mass correlation function between halos of mass 
Ml and M2, separated by a co- moving distance R is (see 
Scannapieco & Barkana 2003 and references therein) 
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W is the window function (top- hat in real space), P{k) the 
power spectrum and p„i is the cosmic mass density. The 
dark-matter halo correlation function for halos of mass M 
is obtained from the product of the mass correlation func- 
tion ^m(M, M, R) and the square of the ratio between the 
variances of the halo and mass distributions. This ratio, 6, 
is defined as the halo bias; its value for a halo mass M may 
be approximated using the Press-Schechter formalism (Mo 
& White 1996), modified to include non-spherical collapse 
(Sheth, Mo & Tormen 2001) 
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where ^c,o (~ 1.69) is the critical overdensity threshold for 

spherical collapse, Sc{z) = 5cfi/D{z), D{z) is the growth 
factor at redshift z, a is the variance on a mass-scale M, 
V = 6l{z)/a'^{M), v' = Vazy, a = 0.707, b = 0.5 and 
c = 0.6. This expression yields an accurate approxima- 
tion to the halo bias determined from N-body simulations 
(Sheth, Mo & Tormen 2001). 

We need the halo correlation function between halos of 
different masses. The bias (equation 3) yields the excess 
probability relative to the mass distribution of finding a 
halo at any point. We may therefore determine the halo 
correlation function for halos of mass Mi and M2 using 
the product of the biases for individual masses, 

eh(Mi, M2, R) = b{Mi, z)b{M2, z)U(Mi, M2, R)D{z)^. 

(4) 

If the correlation function between halos forming at dif- 
ferent redshifts is required, then the more general ana- 
lytic treatment of Scannapieco & Barkana (2003) should 
be used. 

We would like to constrain the relationship between 
SMBH mass, Mbh, and galactic halo mass, Mhaio, by 
constructing a theoretical quasar correlation function for 



comparison with observational data. Observations of 
the quasar correlation function measure the clustering of 
quasars of a certain luminosity, rather than the cluster- 
ing of halos or SMBHs of a particular mass. We therefore 
need to associate the _B-band luminosity Lb of a quasar 
with the halo mass of its host galaxy Mhaio- We begin by 
assuming quasars to have a spectral energy distribution 
corresponding to the median of their population (Elvis et 
al. 1994). We then allow the SMBH to shine at a fraction 
T] of its Edclington luminosity, so that in solar units, the 
B-band luminosity of the quasar is 



-^ = 5.73x10^2 
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We also need to speciiy a relation between SMBH and 
halo mass [Mbh = /(Mhaio), see § 3]. Given a lumi- 
nosity Lb, quasar pairs shining with Eddington fractions 
771 and 772 have their correlation function specified by 
equation (4) where Mi_2 = (Mbh, 1,2) and Mbh, 1,2 = 
Lb/(5.73x lO3 77i,2iB,0)M0. 

The correlation function for a population of quasars with 
a luminosity Lb can therefore be computed by drawing 
pairs of 77- values from an appropriate probability distribu- 
tion, dPohs/drj. The correlation function of quasars with 
luminosity Lb is 



Cq(-^B) = (&(Lb, ??!,%)), 



(6) 



where 771^2 axe drawn from dPohs/d-q and angular brack- 
ets denote an average over the probability distribution of 
?7- values. 

We consider two distributions for rj in this paper. In 
the first case all quasars shine at their Eddington lu- 
minosity, so that dPobs/dr] cx S{ri — 1), where 6 is the 
Dirac delta function. This corresponds to an idealized 
quasar lightcurve that is a tophat-function in time, r]{t) = 
9(— t)8(t — iq), whereby the quasar shines at the Edding- 
ton luminosity (77 = 1) throughout its lifetime, tq. In the 
second case, the observed distribution is assumed be flat in 
the logarithm of 77 so that dPohs/d{logr]) is constant. This 
is the form of distribution expected for a quasar lightcurve 
of the exponential form Q{—t)ri{t) = exp(— t/tq). 

3. THE Mbh-Mhalo RELATION 

Dormant SMBHs are ubiquitous in local galaxies 

(Magorrian et al. 1998) . The masses of these SMBHs scale 
with physical properties of their hosts (e.g. Magorrian et 
al., 1998; Mcrritt & Ferrarese 2001; Tremainc ct al. 2002). 
Motivated by local observations (Ferrarese 2002), we as- 
sume a relation Mbh = / (Mhaio, z) between SMBH mass 
and the mass of the host dark matter halo. This rela- 
tion may be applied to the calculations of the correlation 
function of quasars because the masses of SMBHs scale 
the same way with the physical properties of their host 
galaxies in both quiescent and active galaxies (McLure & 
Dunlop 2002). The function / is constrained by local ob- 
servations (Ferrarese 2002), and we consider two different 
forms for its redshift dependence: 

• Case A: we assume that SMBH mass is correlated 
with the halo circular velocity. This scenario is sup- 
ported empirically by Shields et al. (2003) who stud- 
ied quasars out to 2; 3 and demonstrated that the 
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relation between Mbh and the stellar velocity dis- 
persion does not evolve with rcdshift. This is ex- 
pected if the mass of the black-hole is determined 
by the depth of the gravitational potential well in 
which it resides, as would be the case if growth is 
regulated by feedback from quasar outflows (e.g. 
Silk & Rees 1998; Wyithe & Loeb 2003). Express- 
ing the halo circular velocity, Vc, in terms of the 
halo mass, M, and redshift, z, the redshift depen- 
dent relation between the SMBH and halo masses 
may be written as 



Mbh(Mhaio, z) = const X 
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where e is a constant, ({z) is close to unity and 
defined as C = [(f^™/0^)(A,/187r2)], = [1 + 
{nA/nm){l + z)-^]-\ Ac = 187r2+82rf-39rf2, and 
d = — 1 (see equations 22-25 in Barkana & Loeb 
2001 for more details). 

Ccise B: we consider the scenario in which the 
SMBH mass maintains the same dependence on its 
host halo mass at all redshifts, namely 
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In models All, AIII, BII and Bill, the range 0.01 < 
< 1 is considered. We refer to the above labels when 
presenting our results through the remainder of the paper. 

The models described above utilize a Afbh-Mhaio rela- 
tion with a normalization derived under the simplifying 
assumption of a flat rotation curve. However, the circu- 
lar velocity at the radii probed by observations of galaxies 
is likely to be larger than the circular velocity at their 
virial radius. The normalization e in equations (7) and 
(8) may therefore be larger than the value esis considered 
here (Ferrarese 2002). For any universal halo profile, this 
difference is degenerate with a simple renormalization of 
7]. The results in the following section may be applied to 
a case where the normalization is e' = Ce by substituting 
a correspondingly lower value of the Eddington luminosity 
T]' =rj/C. The allowed range of the product er] is explored 
in § 5. 

4. COMPARISON TO OBSERVATIONS 

Next we compare model correlation functions to obser- 
vations from the 2dF quasar redshift survey (Croom et 
al. 200 lb, c). Figure 1 shows the correlation function ob- 
tained by assuming that SMBH mass scales with halo cir- 
cular velocity and that all quasars shine at their Eddington 
luminosity (case AI). We show the predicted correlation 
function at various redshifts, and several apparent B-band 
magnitudes^, as determined from the relation 



B = 5.5-2.51ogio 



The normalizing constant in these relations has an ob- 
served"* value of e = esiS ~ 10^^'^ at z = (Fer- 
rarese 2002); where SIS denotes the underlying assmnption 
that the halo mass profile resembles a Singular Isothermal 
Sphere. In the following section we construct model cor- 
relation functions assuming both of these scenarios, and 
compare the results to observations. We will show that 
case A is consistent with the data of the 2dF quasar red- 
shift survey, while case B is not. 

3.1. Summary of models 

Before proceeding to the comparison with the observed 
correlation length we label six models of interest: 

• (AI) The fiducial model: Mbh on v^, r] = 1, 

e = esis- 

• (All) Mbh cx vl, dPohs/d{logr]) = const, e = esis- 



- 2.5(1 + a^) logio (1 + 



(9) 



where = —0.5 is the typical slope of the quasar power- 
law continuum (with a frequency dependence of the spec- 
tral flux /jy cx z/"'^), and Dl is the luminosity distance. 
The magnitude limit of the 2dF quasar redshift survey is 
B ~ 20.85, and so the data should be compared with the 
model correlation function corresponding to B ^ 20 — 21 
(solid and dashed lines) . The lower right panel of Figure 1 
shows the 2dF correlation function generated from the en- 
tire quasar sample. This data is compared with a model 
correlation function computed at the median quasar red- 
shift of the survey ((z) ^ 1.5). The B = 20.85 prediction 
at this redshift is replotted in the other panels of figure 1 
to guide the eye (thick gray curves). We find that the 
correlation function depends more sensitively on redshift 
than on luminosity over the ranges under consideration. 

4.1. The correlation function for models with Mbh oc 

Variation of clustering with redshift and apparent mag- 
nitude is demonstrated more clearly in Figure 2. The up- 
per panels show the dependences of the clustering length, 
Ro [defined through the condition ^q(J?o) = 1] for the fidu- 
cial case where the SMBH mass scales with halo circular 
velocity and quasars shine at their Eddington limit (case 
AI). The correlation length is plotted for various values of 
apparent magnitude as a function of redshift (left panel), 
and various redshifts as a function of apparent magnitude 
(right panel). The data for the correlation length as a 

^ We have used the normalization derived by Ferrarese (2002) under the simplifying assumption that the virial velocity of the halo represents 
its circulaj: velocity. 

* We refer to the Johnson B-magnitude, which for a typical quasar is related to the photographic b-band used in the 2dF quasar redshift survey 
by B « & + 0.06 (Goldschmidt & Miller 1998) 
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Fig. 1. — Predicted correlation function at various redshifts, in coinparison to the 2dF data (Croom et al. 2001b). The dark hues show 
the correlation function predictions for quasars of various apparent B-band magnitudes (B). The 2dF limit is B ~ 20.85. The lower right 
panel shows data from entire 2dF sample in comparison to the theoretical prediction at the mean quasar redshift of {z) = 1.5. The B = 20.85 
prediction at this redshift is also shown by thick gray lines in the other panels to guide the eye. This case assumes M^h oc , normalized to 
local observations, and rj = 1 (case AI). 



function of redshift was computed from a sample having 
B < 20.85, drawn from the 2dF quasar redshift survey 
(Croom et al. 2001b, c). The gray error bars show the val- 
ues for Ro in different redshift bins, while the correlation 
length measured for the whole sample (plotted at z = 1.5 
and B = 20.3 in the left and right hand panels respec- 
tively) is shown by the dark error bar. We find the predic- 
tions for Rq to be consistent with observations. In partic- 
ular, both the value of the clustering length and its slow 
evolution with redshift and apparent magnitude are consis- 
tent with current data. There is little variation of the clus- 
tering length with apparent magnitude over the range con- 
sidered, although there is a tendency for brighter quasars 
to be more highly clustered, as was noted by Croom et 



al. (2002). 

Note that the K-correction assumed in equation (9) is 
only valid out to z ~ 3 where Lya absorption enters into 
the -B-band. As a result the extrapolations of redshift de- 
pendence to 2 > 3 shown in Figures 2 and 3 do not include 
Lya absorption in the calculation of the quasar luminos- 
ity. However, these extrapolations arc included to illus- 
trate the evolution of the clustering length at yet higher 
redshifts (as measured in bands redward of the redshifted 
Lya line). 

The lower panels of Figure 2 show the corresponding de- 
pendence of the clustering length on redshift and apparent 
magnitude for quasars that shine at Eddington fractions r] 
lower than unity. Only the cases of B = 20.3, correspond- 
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Fig. 2. — The correlation length, Rq [defined through ^q(i?o) = 1] plotted as a function of redshift (left hand panel) and apparent magnitude 
(right hand panel). The data is from the 2dF quasar redshift survey (Croom et al. 2001b, c), with the gray error bars showing the correlation 
length from different sub-samples, and the dark error bar at z = 1.5 (left hand panels) and B = 20.3 (right hand panels) showing the 
correlation length for the full sample. The theoretical model assumes M^i^ oc v'^. In the upper panels rj = I (case AI), and the curves are 
shown for different apparent magnitudes (left), and different redshifts (right). In the lower panels we show curves corresponding to diflferent 
choices for 17 and e, (cases All and AIII) assuming B = 20.3 and z = 1.5 respectively. 



ing approximately to the median quasar B-magnitude (left 

hand panel) and z = 1.5 (right hand panel) arc shown, and 
the curve corresponding to 77 = 1 (case AI) is reproduced 
from the upper panels for comparison (solid lines). As- 
suming the local normalization e = esis and that rj ranges 
between 0.01 and 1 with a constant rfPobs/'^(log?7) (ease 
All), the typical quasar of apparent magnitude B must 
be powered by a more massive SMBH and hence reside 
in a more massive galaxy than in the rj = 1 case. This 
results in a larger clustering length for quasars of a given 
apparent magnitude B. The resulting evolution, indicated 
by the dashed lines, appears to be incompatible with the 
data. The increase in Rq may be compensated for by an 
increase in the normalization of the Mbh — -^haio rela- 
tion. The dot-dashed curves show the result of increasing 
this normalization relative to local observations by a fac- 
tor of 10 (f = lOesis)- This factor equals the inverse of the 
median of the distribution dP/d{\ogr]) = const. The re- 
sulting curve (corresponding to case AIII) is very similar 
to our fiducial case (77 = 1). 

These results demonstrate that while a spread in the 
values of rj for different quasars results in the correlation 
function being measured over a wider range of halo masses, 
the correlation length is primarily sensitive to the charac- 
teristic halo mass, as specified by the characteristic value 
of 77. Furthermore, the results demonstrate the important 
point that the correlation length measures the product of 



the normalization in the Mbh-Mhaio relation and the me- 
dian fraction of the Eddington luminosity at which quasars 
shine, ?7med- With respect to the correlation function, this 
quantity er^med (which is the same in cases AI and AIII) 
is more fundamental to the correlation function than the 
quasar lifetime since it predicts the correlation length in- 
dependent from consideration of quasar number counts. 
It is important to note that since the normalization of the 
local Mbh — M^aio relation is estimated (e = egis, Fer- 
rarese 2002), the value of the quasar correlation length 
at low redshift implies that quasars spend most of their 
bright phase shining near their Eddington luminosity. 

4.2. The correlation function for models with 

We next examine whether the correlation function ob- 
served by the 2dF survey may be used to constrain the 
evolution of the Mbh-Mhaio relation. In the previous sub- 
section, we have considered the case of a SMBH mass that 
scales with the fifth power of circular velocity independent 
of redshift (equation 7, case A). The alternative case B 
example postulates that the SMBH mass has a redshift 
independent scaling with the halo mass rather than halo 
circular velocity (equation 8). The upper panels of Fig- 
ure 3 plot the correlation length in this case with rj = 1 and 
e = esis (case BI). Both the evolution of the correlation 
length with redshift for different apparent magnitudes (left 
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Fig. 3. — The correlation length, as a function of redshift (left hand panel) and apparent B-band magnitude (right hand panel). The 
data is from the 2dF QRS (Croom eta al. 2001b, c), with the gray error bars showing the correlation length from different sub-samples, and 
the dark error bars at z = 1.5 (left hand panels) and B = 20.3 (right hand panels) showing the correlation length for the full sample. The 
5/3 

model assumes M^h ^^iiAo' upper panels r] = \ (case BI), and the curves are shown for different apparent magnitudes (left), and 

different redsliifts (right). In the lower panels we show curves corresponding to different choices for r) and e (cases BII and Bill), assuming 
S = 20.3 and z = 1.5 respectively. 



hand panel), and the variation of correlation length with a 
given apparent magnitude at various redshifts (right hand 
panel) are shown. The case B Mbh — Afhaio relation results 
in a correlation length that varies significantly more with 
redshift than the data requires. The lower panels show 
examples with values of ry that differ from unity. Only the 
cases of B = 20.3 (left hand panel) and z = 1.5 (right hand 
panel) are shown, and the result assuming t] = 1 (case BI) 
is reproduced from the upper panels for comparison (solid 
lines). If e = esiSj smd rj ranges between 0.01 and 1, with 
dPobs/dr] oc const (case BII), then the typical quasar of 
apparent magnitude B must have a larger clustering length 
(dashed lines). This may be compensated for by increasing 
the normalization of the MBH~-Mhaio relation by a factor 
of 10 (e = lOesis, case Bill) as considered before (dot- 
dashed lines). 

The rapid evolution with redshift that results from a 

case B Mbh Afhaio relation could be diminished if the lu- 
minosities of quasars were a higher fraction of their Ed- 
dington limit at higher redshift. However we have seen 
that the locally observed normalization of the Mbh-Afhaio 
relation requires most present-day quasars to be shining 
near their Eddington luminosity. Quasars at high redshift 
would therefore have to be super-Eddington, shining at 
^ rj ^ (1 + z)^/^ (77 ~ 30 at z ^ 3) in order to repro- 
duce the observed slow evolution of the clustering length. 



This requirement on 77 may be somewhat relaxed in models 
having e > esiS) as is discussed in the next section. 

5. CONSTRAINING T] AND e USING THE QUASAR 
CORRELATION FUNCTION 

The Mbh-Afhaio relation for local galaxies requires a 
choice for the relationship between the circular velocity 
at the radius probed by observations of galaxies, and the 
circular velocity at the virial radius. This relationship de- 
pends on the dark-matter mass profile and is still uncertain 
due to the gravitational interaction between the baryons 
(which cool and often dominate gravity near the center 
of galaxies) and the dark matter. Ferrarese (2002) pre- 
sented Mbh^-Mhaio relations for three different cases, which 
may be summarized in the present notation by e = esiSi 
e = cnfw = 3.7esis, and e = eso2 — 25esis- The second 
case corresponds to the universal profile of Navarro, Frenk 
& White (1996, NEW), and the third to the normalization 
derived based on the weak lensing study of Seljak (2002). 

The results in the preceding sections have been pre- 
sented for fiducial models with (e, r?) = (esiS) 1)) (esiS)O.l) 
and (lOesis, 0.1). It was demonstrated that the correlation 
function is primarily dependent on the product between e 
and the median value of 77, namely er^med- This raises the 
possibility that eTjmed could be constrained using correla- 
tion function data. We have computed the evolution of the 
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Fig. 4.— The probability P(x^) of case A and case B models as a function of feVmed- The dark lines correspond to fits to the redshift 
dependent correlation length Ro{z) measured by the 2dF quasar redshift survey (Croom et al. 2001b). The gray lines show results for a 
hypothetical larger future sample whose errors on Rq are reduced by a factor of \/2 relative the the 2dF quasar survey. 



correlation length Ro{z) for a set of models having a range 
of values for the product er/med = /eesis??med, and calcu- 
lated the reduced x^-statistics as well as the confidence 
with which a model is excluded by the 2dF data, P(x^). 
In figure 4 we plot the resulting P{x^) as a function of 
feVmed (dark lines). 

The results suggest that Case A models (solid lines) fa- 
vor /eJymed ^ 1^5 indicating that e is a factor of a few 
larger than egis for ry ^ 1, as would be expected from the 
universal NFW profile, or that rymed has a maximum value 
of a few tenths. If e = eso2, case A requires r?med ~ 0.08. 
Figure 4 shows that Case B models (dashed lines) are more 
restricted, requiring /erymod ^ 10 — 30. This implies that 
e = esis and e = eNFW cannot be accommodated since 
the required values of rjmed ^ 10 and 77med ^ 3 respec- 
tively, are considered unphysical in most accretion models 
(Nityananda & Narayan 1982; see, however, Ruszkowski & 
Begelman 2003). The more extreme case of e = eso2 may 
be marginally accommodated by case B; however even in 
this case the quasars must shine somewhat in excess of 
their Eddington limit. 

Figure 4 also shows P(x^) as a function of /e?ymcd for a 
hypothetical future survey that would yield the same evo- 
lution of -Ro with redshift, but with errors that are smaller 
by a factor of \/2 (light lines). Future surveys will better 
constrain the gradient dRo/dz, and help discriminate fur- 
ther between different models of the Afbh Afhaio relation. 

Our discussion so far has been independent of the quasar 
luminosity function. In the next section we discuss how the 
luminosity function may be used to determine the redshift 
evolution of the quasar lifetime. 

6. THE LUMINOSITY FUNCTION AND QUASAR LIFETIME 

In the previous sections we have demonstrated that the 
evolution of the Mbh^Afhaio relation with redshift may be 
constrained directly from the quasar correlation function. 
We have shown that the critical parameter determining 
the SMBH mass is the halo circular velocity rather than 



the halo mass. We reiterate that this result does not 
rely on the luminosity function of quasars. On the other 
hand, the number counts of quasars are proportional to 
the quasar lifetime. Therefore, the quasar lifetime and 
its evolution with redshift may be determined by com- 
paring the observed quasar number counts to the SMBH 
mass function. The latter may be obtained from the Press- 
Schechter (1976) mass-function of dark-matter halos (with 
the improvement of Sheth & Tormen 1999) combined with 
the Mbh-Mhaio relation. 

Theoretical models of the quasar luminosity function 
often associate quasar activity with major mergers (e.g. 
Kauffmann & Haehnelt 2000; Wyithe' & Loeb 2003). 
A simple model that attributes quasar activity to ma- 
jor mergers, assumes a Mth-Mhaio relation that is de- 
scribed by equation (7), and sets the quasar lifetime equal 
to the dynamical time of a galactic disk, accurately re- 
produces the entire optical and X-ray luminosity func- 
tions of quasars at redshifts between 2 and 6 (Wyithe & 
Loeb 2003). At lower redshifts (z < 2), this prescription 
correctly predicts the luminosity function at luminosities 
below the characteristic break (Boyle et al. 2000). In dif- 
ference to Wyithe & Loeb (2003) we have assumed the 
quasar activity to take place in the larger SMBH before 
coalescence of the SMBH binary formed during the merger 
(Volonteri et al. 2002). Figure 5 shows the comparison of 
this model with observations. Specifically, we show the 
fiducial case Al (which corresponds to the model pre- 
sented in Wyithe & Loeb 2003) using equation (7) with 
e = esis and r] = \ (solid curves). We have shown that 
case Al describes the observed quasar correlation function 
(see figures 1-2). However as discussed in § 5, case-B mod- 
els for the Mbh-Mhaio relation also provide acceptable fits 
to R(}{z) so long as e > lOesis and ?7 ~ 1. Therefore a 
corresponding luminosity function model was also gener- 
ated for this case assuming a constant lifetime of = 10^ 
years. We have plotted this model in figure 5 (dot-dashed 
lines). 
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Fig. 5. — Compajrison of the model and observed quascir luminosity functions at various redshifts. The data points at z < 4 are summarized 
in Pei (1995), while the light lines show the double power-law fit to the 2dF quasar luminosity function (Boyle et al. 2000). At z ~ 4.3 and 
z ~ 6.0 the data is from Fan et al. (2001a;2001b;2003). The gray regions show the 1-cr range of logarithmic slope ([—2.25, —2.75] at z ~ 4.3 
and [—1.6, —3.1] at z ~ 6), and the vertical bars show the uncertainty in the normalization. The open circles show data points converted from 
the X-ray luminosity function (Baxger et al. 2003) of low luminosity quasars using the median quasar spectral energy distribution (Elvis et 
al. 1994). The model luminosity functions shown correspond to cases Al (solid curves) and a model assuming case B evolution with e = eso2, 
ri = 1 and a constant lifetime of tq = 10^ years (dot-dashed curves). The dashed pairs of vertical lines show the range corresponding to 
apparent B-magnitudes between 18.25 and 20.85. 



We find that our fiducial model (case Al; solid lines) 
with a quasar lifetime that equals the dynamical time of a 
galactic disk, reproduces the observed luminosity function 
over a wide range of redshifts and luminosities (see Wyithc 
& Loeb 2003). Our second model (case B, e > lOesiS) 

= 1, = 10''; dot-dashed lines) also reproduces the ob- 
served luminosity function at z < 4.5. However at ~ 6 
this model underestimates the observed number counts by 
two orders of magnitude. Case B evolution for the Mbh^ 
Mhaio relation is therefore challenged by the existence of 
the highest redshift quasars, which would imply a quasar 
lifetime that exceeds the age of the universe at z ~ 6. 
Moreover, given the Salpeter (1964) time of ^ 4 x 10^ 
years, the SMBHs powering these quasars would grow to 
25 times their initial mass during their lifetime. Provided 
that quasar activity is trigged by major mergers, the re- 
quirement that the Mbh Mhaio relation reproduce the ob- 
served evolution in quasar clustering, together with the 
observed evolution in the number counts of quasars there- 
fore implies case A evolution, and a quasar lifetime that is 
set by the dynamical time of the host galaxy. We reinforce 
the important qualitative result that since the evolution of 



quasar correlations require SMBHs to take a larger fraction 
of the halo mass at high redshift (case A, see equation 7), 
the observed evolution in the number counts requires a 
quasar lifetime that is shorter at higher redshifts, scaling 
approximately as the dynamical time of a dark matter halo 
[oc il + z)-^/% 

We point out that the model luminosity functions fail 
to reproduce the rapid drop in the density of luminous 
quasars seen at low redshift. This drop is thought to be 
due to the consumption of cold gas in galaxies by star 
formation (e.g. Kauffmann & Haehnelt 2000) and the in- 
hibition of accretion onto massive halos at late times (e.g. 
Scannapieco & Oh 2004). However, most of the quasars 
in the 2dF quasar redshift siu'vcy have i3-magnitudes that 
lie below the characteristic break in the luminosity func- 
tion at low redshift (Boyle et al. 2000; Groom et al. 2001). 
This can be seen in Figure 5 where the vertical dashed 
lines show the luminosities bracketing the apparent mag- 
nitude range 18.25 < B < 20.85. Since the model correctly 
predicts the luminosity function within most of this lumi- 
nosity range at all redshifts, we are justified in deriving 
combined constraints from the luminosity and correlation 
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functions of quasars. 

7. DO BRIGHT LOW-REDSHIFT QUASARS SHINE NEAR 
THEIR EDDINGTON LUMINOSITY? 

In § 4 we demonstrated that the 2dF quasar sample 

has a correlation length that is consistent with the local 
Mbh — Afhaio relation, and that quasars shine near their 
Eddington luminosity. In addition it was shown that this 
consistency extends to both bright and faint subsamples of 
quasars. The brightest 2dF quasars (18.1 < B < 19.9) pre- 
sented by Croom ct al. (2002) arc mostly brighter than the 
characteristic break luminosity in the differential quasar 
number comits (e.g. Boyle et al. 2000). These bright 
quasars have a lower abundance than expected from merg- 
ing halos at z < 2 (sec figure 5). Either the lifetime (or 
occupation fraction) of bright quasars is lower than ex- 
pected at z < 2 or else these quasars shine well below 
their Eddington limit. Figure 2 indicates that the slow 
dependence of Ro with B reported by Croom et al. (2002) 
is consistent with the expected trend if all quasars shine 
near their Eddington luminosity. This consistency sug- 
gests that bright quasars are under represented at z < 2 
because of a reduced quasar lifetime (or a reduced occu- 
pation fraction) rather than a reduced value of r] during 
their brightest phase. 

8. DISCUSSION 

In this paper we have demonstrated that the relation- 
ship between the mass of a SMBH and its host dark-matter 
halo is the fundamental quantity that determines the clus- 
tering of quasars. In particular the Mbh-Mhaio relation 
and its evolution with redshift, may be determined directly 
from the evolution in the correlation length of quasars, in- 
dependent of consideration of the quasar luminosity func- 
tion or assumptions about the quasar lifetime. Beginning 
with the locally observed Mbh Mhaio relation, we have 
shown that the observc;(l correlation length of quasars in 
the 2dF quasar redshift survey is consistent with SMBHs 
that shine at ~ 10-100% of their Eddington luminosity 
during their bright quasar episode. Moreover, the evo- 
lution of the clustering length with redshift is consistent 
with a SMBH mass whose redshift-independent scaling is 
with the circular velocity of the host dark matter halo 
(see equation 7). In contrast, it appears that a relation 
for the SMBH mass whose redshift-independent scaling is 
with the mass of the host dark matter halo (equation 8) 
would have resulted in too much evolution of the clustering 
length with redshift. 

The portion of the 2dF quasar survey from which the 
data used in this study were drawn, contains ^ 10^ 
quasars. Upon completion, the Sloan Digital Sky Sur- 
vey will have spectra for ten times this number of quasars, 
spread over a larger redshift range (York et al. 2000). This 
will allow more accurate determination of the correlation 
length and its evolution, and hence provide more accu- 
rate constraints the SMBH population. The selection of 
quasars for the spectroscopic survey of SDSS at 2; < 3 is 
restricted to ?'* < 19. Thus, most of the low redshift SDSS 
quasars will fall at luminosities that are brighter than the 
characteristic break in the luminosity function (e.g. Boyle 
et al. 2000). These bright quasars have a lower abundance 
than expected from the number of halos merging at 2: < 2. 



The reason for this low abundance may be either that the 
lifetime (or occupation fraction) of bright quasars is lower 
than expected at 2; < 2, or alternatively that the fraction 
of the Eddington limit at which the quasars shine is small. 
The 2dF quasar sample straddles the luminosity function 
break. We have shown that the similarity in the cluster- 
ing statistics of sub-samples of quasars having different 
ranges in luminosity suggest that the lifetime (or occupa- 
tion fraction) of bright quasars is lower than expected at 
z < 2, but that these rjuasars also shine near their Edding- 
ton luminosity. Measurements of the correlation length as 
a function of quasar luminosity in SDSS can help further 
distinguish between these two possibilities. 

As we have shown in a previous paper (Wyithe & 
Loeb 2003), feedback-regulated growth of SMBHs during 
active quasar phases implies Mbh oc v^, consistent with 
the inferred relation for galactic halos in the local universe 
(Ferrarese 2002). Assuming ry = 1 as suggested by obser- 
vations of both low and high- redshift quasars (Floyd 2003; 
Willott, McLure & Jarvis 2003), we have shown that only 
~ 7% of the Eddington luminosity output needs to be 
deposited into the surrounding gas in order to unbind it 
from the host galaxy over the dynamical time of the sur- 
rounding galactic disk. The power-law index of 5 in the 
M\,h-Vc relation, is larger than the values of 4-4.5 inferred 
from the local relation between Mbh and the stellar ve- 
locity dispersion cr^, (Merritt & Ferrarese 2001; Tremaine 
et al. 2002); the difference originating from the observa- 
tion that the Vc-Ci, relation is shallower than linear for 
stellar bulges embedded in cold dark matter halos (Fer- 
rarese 2002). Feedback regulated growth also implies that 
the Mbh-'^c relation should be independent of redshift, as 
implied by the slow evolution in the clustering length of 
quasars. 

A very simple model can therefore be constructed to de- 
scribe the correlation length of quasars and its evolution 
in the 2dF survey, as well as the number counts of quasars 
out to very high redshift {z ~ 6). The model includes 
four simple assumptions, but no free parameters: (i) the 
locally observed Mbh-Mhaio relation extends to high red- 
shifts through equation (7) with the SMBH mass scaling as 
halo circular velocity to the 5-th power; (ii) SMBHs shine 
near their Eddington luminosity with a universal spectrum 
(Elvis et al. 1994) during luminous quasar episodes; (iii) 
quasar episodes are associated with major galaxy mergers; 
and (iv) the quasar lifetime is set by the dynamical time 
of the host galaxy. 

The assumption of a SMBH mass that scales with halo 
circular velocity independently of redshift is supported by 
the observations of Shields et al. (2003) that there is no 
evolution in the M^h-cr* relation out to z ~ 3. The as- 
sumption that quasars shine near their Eddington limit is 
supported by observations of high and low redshift quasars 
(Floyd 2003; Willott, McLure & Jarvis 2003). Overall, the 
sample of available data on the clustering properties and 
the number counts of quasars is most readily explained by 
a quasar lifetime that is set by the dynamical time of the 
host galaxy rather than by the Salpeter (1964) e-folding 
time for the growth of its mass. 

This work was supported in part by NASA grant NAG 
5-13292, and by NSF grants AST-0071019, AST-0204514 
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